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Within the framework of an isospin-dependent quantum molecular dynamics model, 
the rotational flow in reactions of 58 Fe + 58 Fe and 58 Ni + 58 Ni at 40 MeV/nucleon for 
different impact parameters is investigated by analyzing the mid-rapidity azimuthal 
distribution. The rotational observables are also calculated semiquantitatively. For 
the first time, it is found that the more neutron-rich system ( 58 Fe + 58 Fe) exhibits 
stronger rotational collective flow. This isospin dependence of rotational collective flow 
is more appreciable in semi-peripheral collisions and it is shown to mainly result from 
the isospin dependence of nucleon-nucleon cross section rather than the symmetry 
energy. Meanwhile, it is indicated that the rotational flow depends strongly on the 
impact parameter. 

PACS number(s): 25.70.-z, 25.75.Ld, 02.70.Ns 

The azimuthal distribution in heavy ion collisions (HIC's) has proved very useful in the study of 
nuclear equation of state (EOS) as well as reaction dynamics because of its sensitivity to collective 
flow pj]. In particular, the rotational flow may still exist even though the directed flow has disappeared 
|^| . With the recent advance in radioactive nuclear beam physics, one can investigate isospin degrees 
of freedom in nuclear reactions at wide energy ranges for different projectile-target combinations (^^]. 
Therefore, it is very significant to explore the isospin effects on the azimuthal distribution through 
studying HIC's induced by the radioactive nuclei since it may provide people a better probe of isospin 
dependent reaction dynamics. 

We report here results of the first theoretical study on the isospin dependence of the azimuthal 
asymmetry rotational flow from reactions of 58 Fe + 58 Fe and 58 Ni + 58 Ni at 40 MeV/nucleon for 
different impact parameters within the framework of an isospin dependent quantum molecular dy- 
namics (IQMD) model which includes the symmetry energy, Coulomb interaction, isospin-dependent 
experimental N-N cross sections, and particularly the isospin dependent Pauli blocking [||||. In the 
present calculations, the so-called soft EOS with an incompressibility of K = 200 MeV is used and 
the symmetry strength C — 32 MeV without particular consideration || . In the initialization process 
of the IQMD model, the neutron and proton are distinguished from each other Q and meanwhile the 
nonphysical rotations in the initialized nuclei have been removed. 

In the IQMD model simulation, the reaction plane is known a priori and it is defined as the x-z 
plane (z-axis corresponds to the beam direction). The azimuthal angle with respect to the reaction 
plane can be written as 

<f> = arctan(P y /P x ). (1) 

In the present calculations, the calculated results have adopted the average values from 200 to 300 
fm/c to accumulate the numerical statistics since the azimuthal distributions have been well stable 
after 200 fm/c and 200 events are simulated for each impact parameter. Fig. 1 displays the IQMD 
model predicted normalized azimuthal distribution dN/dcf) for all nucleons from 58 Fe + 58 Fe (solid 
circles) and 58 Ni + 58 Ni (open circles) at 40 MeV/nucleon and impact parameter 6 = 6 fm. Figures 
1 (a), (b) and (c) correspond to the target-like rapidity region, mid-rapidity region and projectile- 
like rapidity region, namely, the values of reduced center-of-mass (cm.) rapidity (y/y pr oj)c.m. of 
nucleons considered belong to [—1.5,-0.5], [—0.5,0.5] and [0.5,1.5], respectively. The errors shown 
are statistical. Also included in Fig. 1 are the Legendre polynomial fits up to the second order for the 
resulting azimuthal distribution, i.e., 



dN/d(f) = c(l + a\ cos{4>) + a 2 cos(20)). 



(2) 



The solid (dashed) line is the result of fit for 58 Fe + 58 Fe ( 58 Ni + 58 Ni). In Eq. (2) the coefficient <2i 
represents the strength of the directed flow with preferred emission at <f> — ±180° (i.e., ai < 0) for low 
rapidity values (in the backward hemisphere) and = 0° (i.e., a\ > 0) for high rapidity values (in the 
forward hemisphere) if the incident energy is above the balance energy. The coefficient a 2 reflects a 
flattening of the ellipsoid. A negative value of a 2 (i.e., the azimuthal distribution peaks at <f> = ±90°, 
simultaneously) reflects the squeeze-out effects and a positive value (i.e., the azimuthal distribution 
peaks at <fi — 0° and ±180°, simultaneously) the rotational collective motion. One can see from Fig. 
1 that for low (target-like) rapidity region the azimuthal distribution only peaks at <j> = 0° while it 
only peaks at ±180° for high (projectile-like) rapidity region, which imply that these reaction systems 
display negative deflections (corresponding to negative scattering angles) since the incident energy is 
below the balance energy 0] . 

It is indicated in Fig. 1 (b) that the azimuthal distribution peaks at <fi = 0° and ±180° simultane- 
ously for mid-rapidity nucleons, which are indicative of a rotation-like behavior. The strength of the 
rotational flow, a 2 , equals to 0.20151 and 0.18047 for 58 Fc + 58 Fe and 58 Ni + 58 Ni, respectively. In 
order to observe the impact parameter dependence of rotational flow and explore its isospin depen- 
dence at different impact parameters, we show values of a 2 for mid-rapidity nucleons in Fig. 2 as a 
function of impact parameter b for 58 Fe + 58 Fe (solid circles) and 58 Ni + 58 Ni (open circles) at 40 
MeV/nucleon. The errors shown are statistical. A strong isospin dependence is observed in Fig. 2, 
namely, the more neutron-rich system 58 Fe + 58 Fe displays stronger collective flow than the system 
n8 Ni + 58 Ni, especially in the semi-peripheral collisions (b — 6 fm). In addition, one can see that 
the rotational flow increases with increment of impact parameter b from near central to peripheral 
collisions and then decreases, and finally disappears in the most peripheral collisions since at such 
large impact parameter the interaction between the projectile and target is very weak. 

In order to see more clearly the rotational flow and its isospin dependence, one can calculate 
some characteristic physical quantities of the rotation for mid-rapidity nucleons, such as the angular 
momentum, the moment of inertia, the angular velocity and its rotational energy, with the help of 
classical mechanics. In the present IQMD model calculations, we simply consider the nucleon system 
as a rigid body at every time point and then the angular momentum about the normal of the reaction 
plane (y-axis) can be calculated as follows, 

L V = E = ^2( Z iPxi - x iPzi), ( 3 ) 

where Xi, Zi, p Xi , p Zi are components on the reaction plane of and Pi in the cm. system. The 
moment of inertia about the y-axis can be written as 

Therefore, the angular velocity around the y-axis can be expressed as 

u y = Ly/Iy, (5) 

and its collective rotational energy as 

Erot = \lyU 2 y. (6) 

Nevertheless, the semiquantative description above may provide us some intuitive information about 
the in-plane collective rotational behavior although it is classical and very rough for HIC's at inter- 
mediate energies. Fig. 3 illustrates the time evolution of the absolute value of angular momentum 
per nucleon \L y /A\, moment of inertia per nucleon I y /A, absolute value of angular velocity \to y \, and 
rotational energy per nucleon E rot /A from 58 Fe + 58 Fe (solid circles) and 58 Ni+ 58 Ni (open circles) at 
40 MeV/nucleon and b = 4 fm for mid-rapidity nucleons. The errors shown are statistical. Meanwhile, 
also included in Fig. 3 (a) are the time evolution of the absolute value of angular momentum per 
nucleon around the y-axis for the whole system of reactions 58 Fe+ 58 Fe (solid squares) and 58 Ni+ 58 Ni 
(open squares) and its amplified window figure for clarity, which indicate that the conservation of 
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angular momentum is destroyed slightly, namely, the angular momentum per nucleon for 58 Fe+ 58 Fe 
(58 Ni+ 58 Ni ) conges slightly, i.e., from 1.388 (1.388) ft/nucleon to 1-348 (1.343) ft/nucleon in thc 
whole reaction process, and meanwhile the isospin dependence and the statistical errors are very 
small. The values for 58 Ni + 58 Ni have been offset with 2 fm/c in the horizontal direction for clarity. 
In fact, both L y /A and uj v are negative by Eqs. (3) and (5) since the reactions display negative deflec- 
tions. From Fig. 3 (a), one can see clearly that the increment of angular momentum becomes more 
slowly after about 125 fm/c and tends to its asymptotic value after a longer time when the reaction 
system may have not any further interaction. The moment of inertia reaches its minimum at about 
40 fm/c when the reaction systems reach their maximum overlap and then increases with the reaction 
systems expanding. The angular velocity peaks at about 50 fm/c when the system has just undergone 
the strongest exchanges of all kinds of degrees of freedom and then becomes saturate after about 200 
fm/c. Similarly, E rot /A has maximum at about 100 fm/c and then becomes saturate after about 200 
fm/c. From above analysis, one can conclude that there actually exists in-plane collective rotational 
motion for mid-rapidity nucleons at such incident energy and impact parameter. The similar results 
are observed for other impact parameters. More importantly here, one can find from Fig. 3 that the 
more neutron-rich system 58 Fe + 58 Fe exhibits clearly stronger rotational energy and angular momen- 
tum than system 58 Ni + 58 Ni, which confirms again the conclusion obtained in Fig. 2. Meanwhile, 
these semiquantative calculations also support the view about extracting the rotational flow from the 
azimuthal distribution analysis. 

The isospin dependence of the rotational flow may be a result of the competition among several 
mechanisms in the reaction dynamics, such as the symmetry energy, isospin dependent JV-JV cross 
sections, Coulomb energy, the surface properties of the colliding nuclei, and so on. Here, we investigate 
the influence of the symmetry energy and isospin dependent JV-JV cross section on the rotational flow. 

Using different symmetry energy strength C and parametrizations of JV-JV cross sections, we show 
in Fig. 4 the IQMD model predicted azimuthal distribution from 58 Fe + 58 Fe (solid circles) and 
58 Ni + 58 Ni (open circles) at 40 MeV/nucleon and b = 4 fm for mid-rapidity nucleons. The errors 
shown are statistical. Meanwhile, the results of Legendre polynomial fits according to Eq. (2) for 
58p e _|_ 58p e (goijd line) and 58 Ni + 58 Ni (dashed line) as well as the resulting a\ and a 2 are also 
included in Fig. 4. In Fig. 4 (a) we use C = (no symmetry energy) and Cugnon's N-N cross section 
<7cug which is isospin independent. It is indicated that a-i of 58 Fe + 58 Fe is slightly less than that 
of 58 Ni + 58 Ni, which may result from the difference of Coulomb interaction and surface properties 
between the two systems. The case of using C — 32 McV and ocug, is plotted in Fig. 4 (b). For the 
results shown in Fig. 4 (c) we use C = 32 MeV and experimental JV-JV cross section a cxp which is 
isospin dependent. One can see from Fig. 4 that the values of a\ are very small since the transverse 
momentum disappears at mid-rapidity and meanwhile can find that the symmetry energy enhances 
the rotational flow (a 2 ) while cr cxp reduces it more strongly. It is also indicated that the influence of 
(T oxp on system 58 Ni + 58 Ni is stronger than that on the system 58 Fe -I- 58 Fe, which just results in 
the observed isospin dependence. This is easy to understand since the neutron-proton cross section 
is about three times larger than the neutron-neutron or proton-proton cross section for a cxp , which 
results in more JV-JV collisions for 58 Ni + 58 Ni. From above analysis, one can conclude that the isospin 
dependence of rotational flow seems to mainly result from the isospin dependence of JV-JV cross section 
and the symmetry potential has little influence on it. 

In summary, by using the IQMD model, we studied for the first time the in-plane rotational flow 
in reactions of 58 Fe + 58 Fe and 58 Ni + 58 Ni at 40 MeV/nucleon for different impact parameters by 
analyzing the mid-rapidity azimuthal distribution and calculating semiquantitatively the rotational 
observables. A strong isospin dependence of the rotational flow has been found, namely, the more 
neutron-rich system exhibits stronger rotational flow. This isospin dependence is more appreciable in 
semi-peripheral collisions and it is shown to mainly result from the isospin dependence of JV-JV cross 
section rather than the symmetry energy. Meanwhile, it is indicated that the rotational flow depends 
strongly on the impact parameter, namely, it increases with increment of impact parameter from 
near central to peripheral collisions and then decreases, and finally disappears in the most peripheral 
collisions. 
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FIGURE CAPTIONS 

Fig. 1 The IQMD model predicted normalized azimuthal distribution dN/dcf) for all nucleons from 
58 Fe + 58 Fe (solid circles) and 58 Ni + 58 Ni (open circles) at 40 MeV/nucleon and impact param- 
eter b — 6 fm for different rapidity regions, i.e., the target-like rapidity region (a), mid-rapidity 
region (b) and projectile- like rapidity region (c). The solid (dashed) line represents the Legendre 
polynomial fit up to second order for the calculated result of 58 Fe + 58 Fe ( 58 Ni + 58 Ni). 

Fig. 2 The values of ai for mid-rapidity nucleons as a function of impact parameter b for 58 Fe + 58 Fe 
(solid circles) and 58 Ni + 58 Ni (open circles) at 40 MeV/nucleon. The lines are plotted to guide 
the eye. 

Fig. 3 The time evolution of the absolute value of angular momentum per nucleon \L y /A\ (a), moment 
of inertia per nucleon I y /A (b), absolute value of angular velocity \w y \ (c), and rotational energy 
per nucleon E rot /A (d) from 58 Fe + 58 Fe (solid circles) and 58 Ni + 58 Ni (open circles) at 
40 MeV/nucleon and b = 4 fm for mid-rapidity nucleons. Also included in (a) are the time 
evolution of the absolute value of angular momentum per nucleon for the whole system of 
reactions 58 Fe+ 58 Fe (solid squares) and 58 Ni+ 58 Ni (open squares) and its amplified window 
figure. 

Fig. 4 The IQMD model predicted azimuthal distribution for 58 Fe + 58 Fe (solid circles) and 58 Ni 
+ 58 Ni (open circles) at 40 MeV/nucleon and b = 4 fm for mid-rapidity nucleons by using 
different symmetry energy strength C and paramctrizations of JV-JV cross sections: C = (no 
symmetry energy) with Cugnon's N-N cross section trcug (a), C — 32 (MeV) with ccug (b), and 
C = 32 (MeV) with experimental N-N cross section a cxp (c). Meanwhile, the results of Legendre 
polynomial fits according to Eq. (2) for 58 Fc + 58 Fe (solid line) and 58 Ni + 58 Ni (dashed line) 
as well as the resulting a\ and a-i are also included. 
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